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Abstract:  We  investigate  the  transient  temperature  rise  in  a  composite  sandwich  structure  induced  by  a 
stationary,  dithering,  or  rotating  laser  beam.  We  restrict  our  study  to  the  composite  sandwich  structure 
with  carbon  fiber  as  skin  materials  and  honeycomb  as  core  materials.  Our  numerical  simulations  indicate 
that  the  maximum  temperature  rise  behaves  as  a  quadratic  function  of  the  reciprocal  of  the  skin  thickness 
or  the  effective  beam  size. 
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1.  Introduction 

In  recent  years  developing  materials,  such  as  materials  of  composite  sandwich 
structure,  to  enhance  the  damage  tolerance  of  the  object  has  become  more  and  more 
important.  The  goal  of  this  paper  is  to  investigate  the  maximum  temperature  rise  of  a 
composite  sandwich  structure  induced  by  a  stationary,  rotating  and  dithering  laser  beam 
and  to  seek  an  optimal  design  of  the  composite  sandwich  structure. 

The  research  on  the  temperature  rise  induced  by  a  laser  beam  is  not  new.  In  the 
past  decades  many  theoretical  and  numerical  studies  of  temperature  profiles  induced  by 
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laser  radiation  in  solids  were  conducted  (Cline  and  Anthony,  1977;  Lax,  1977  and  1978; 
Bertolotti  and  Sibilia,  1981;  Burgener  and  Reedy,  1982;  Calder  and  Sue,  1982;  Moody 
and  Hendel,  1982;  Sanders,  1984;  Araya  and  Gutierrez,  2006;  Sistaninia  et  al.,  2009). 
Temperature  distributions  in  a  two-layer  structure  by  a  scanning  laser  beam  were 
investigated  in  (Burgener  and  Reedy,  1982).  Recently,  detailed  studies  of  the 
temperature  rise  induced  by  a  rotating  or  dithering  laser  beam  were  given  on  a 
semi-infinite  domain  (Zhou,  201  la),  a  solid  with  finite  geometry  (Zhou  and  Tan,  201 1 ; 
Tan  and  Zhou,  2012),  a  two-layer  structure  (Zhou,  2011b;  Zhou,  2012a)  and  a 
composite  with  a  sandwich  structure  (Zhou,  2012b). 

Composite  materials  are  man-made  materials  and  they  are  fabricated  by 
attaching  two  stiff  but  thin  skins  to  a  thick,  lightweight  core.  The  thick  core  provides  the 
sandwich-structured  composite  with  high  bending  stiffness  and  overall  low  density.  In 
(Zhou,  2012b),  the  performance  of  a  sandwich-structured  composite  with  different 
combinations  of  skin  materials  and  core  materials  was  analyzed.  It  was  found  that 
carbon  fibers  are  better  skin  materials  than  E-glass  fibers  and  composite  sandwich 
structures  with  carbon  fiber  skin  and  honeycomb  core  can  serve  extremely  well  as 
heat  shields.  However,  the  scaling  behavior  of  the  maximum  temperature  rise  versus 
different  design  parameters  was  left  unexplored  in  (Zhou,  2012b).  In  this  paper,  we 
extend  our  earlier  preliminary  studies  to  investigate  the  effects  of  design  parameters, 
including  the  thickness  of  the  skin  materials,  the  thickness  of  the  core  materials  and  the 
beam  size,  on  the  maximum  temperature  rise  of  the  composite.  Intuitively,  the  thicker 
the  skin  is,  the  smaller  the  maximum  temperature  rise  will  be  and  it  would  provide  more 
heat  shields.  However,  thicker  skin  will  increase  the  density  of  the  composite  and 
eventually  lose  the  attractive  lightweight  property.  So  there  is  a  compromise  between 
heat  protection  and  lightweight.  We  wish  to  investigate  the  dependence  of  the 
maximum  temperature  rise  of  the  structure  on  the  thickness  of  the  skin  materials  and 
provide  insight  on  the  design  of  the  composite.  Similar  argument  applies  to  the 
thickness  of  the  core  and  we  need  to  investigate  its  effect  as  well.  In  this  paper  we 
limit  our  study  to  a  sandwich-structured  composite  made  of  honeycomb  core 
sandwiched  between  two  thin  layers  of  carbon  fibers  since  this  selection  of  materials 
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has  been  shown  to  provide  an  effective  heat  shield  from  our  previous  studies  (Zhou, 

2012b). 

We  organize  our  paper  as  follows.  In  Sections  2-3  we  present  the  modeling  and 
numerical  solutions  of  the  temperature  distributions  in  a  three-dimensional  composite 
sandwich  structure  induced  by  a  stationary,  dithering  and  rotating  laser  beam, 
respectively.  Finally,  we  draw  conclusions  in  Section  4. 

2.  MATHEMATICAL  FORMULATIONS 

We  first  briefly  review  a  time-dependent  mathematical  model  of  a 
three-dimensional  composite  sandwich  structure  with  one  core  material  sandwiched 
between  two  thin  layers  of  other  materials.  Laser  beams  hit  the  top  surface  of  the 
sandwich  structure.  Figure  1  shows  a  schematic  diagram  of  the  composite  sandwich 
structure.  We  divide  the  finite  structure  into  three  regions:  Region  1  is  a  thin  layer  of 

thickness  dl  in  the  z  direction,  which  is  on  a  substrate  of  dissimilar  material  of 

thickness  d2  (Region  2),  and  Region  3  is  another  thin  layer  of  thickness  di  at  the 
bottom  of  the  structure.  The  materials  in  both  Region  1  and  Region  3  are  the  same. 
Usually  we  set  dx  =  d3  so  the  top  and  bottom  skin  layers  share  the  same  thickness. 


Region  3  (skin) 


Figure  1:  A  schematic  diagram  of  a  three-dimensional  sandwich-structured 
composite. 


580 


HONG  ZHOU,  JOSIAH  PICKETT,  AND  AURELIA  MINUT 


The  temperature  distribution  in  the  composite  sandwich  structure  can  be 
modeled  by  the  standard  non-homogeneous  heat  equation  (John,  1981).  More 
specifically,  in  Region  1  the  governing  equation  is 
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dt 


=  a 


T,  1 
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where  (x,y,z,t)  denotes  the  temperature  rise  at  position  (x,y,z)  and  time  t,  aT1 


is  the  thermal  diffusivity  of  the  material  in  Region  1 ,  KT  l  is  the  thermal  conductivity, 

and  q(x,y,z,t)  is  the  energy  distribution  of  the  incident  laser  beam.  In  the  simulations 

presented  here  we  model  the  distribution  of  a  dithering  or  rotating  Gaussian  beam  as 

q(x,y,z,t)  =  f(x,y,t)8(z-z0).  (2) 


Here  /  has  the  expression 
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where  (xc  (t),yc  (t),z0)  describes  the  position  of  the  dithering  or  rotating  Gaussian 
beam,  (x0,y0,z0)  gives  the  initial  position  of  the  moving  laser  beam,  I0  is  the 
intensity  of  the  laser  beam,  d  is  the  effective  radius  of  the  laser  beam,  and  S(z-z0) 

is  the  Dirac  delta  function.  Equation  (3)  represents  a  rotating  laser  beam  when  a  =  b 
and  a  dithering  beam  when  either  a  or  b  is  zero.  When  a  =  b  =  0 ,  equation  (3) 
simply  represents  a  non-moving  stationary  Gaussian  beam. 

In  Region  2,  where  the  material  has  the  thermal  diffusivity  aT  2  and  thermal 

conductivity  KT  2  ,  the  governing  equation  for  the  temperature  rise  u2(x,y,z,t)  is 
given  by 
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Similarly,  in  Region  3  the  equation  for  the  temperature  rise  u3  (x,y,z,t)  is 
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At  the  interfaces  between  different  regions/materials,  we  impose  the  physical  condition 

that  the  energy  is  conserved  for  the  heat  flow  across  the  interfaces: 
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at  the  interface  between  Regions  1  and  2 
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at  the  interface  between  Regions  2  and  3 


The  initial  condition  for  the  temperature  rises  ux  [x,y,z,t),u2  (x,y,z,t),u3  (x,y,z,t) 

is  zero  to  reflect  the  fact  that  initially,  the  composite  sandwich  structure  shares  the  same 
temperature  as  the  ambient.  The  boundary  conditions  at  the  air/material  interface 
require  that  the  sandwich  structure  is  insulated  at  all  edges.  The  boundary  conditions 
take  the  assumption  that  no  energy  leaks  into  the  ambient  at  the  air/material  interface. 
This  is  fairly  a  good  approximation  for  most  materials  under  consideration  because  heat 
flow  by  conduction  through  the  material  is  usually  much  larger  than  heat  loss  by 
radiation  or  convection  at  the  air/material  interface. 

For  the  top  (Region  1)  and  bottom  (Region  3)  layers  of  the  composite  sandwich 
structure,  the  materials  are  chosen  to  be  carbon  fibers.  Carbon  fibers  are  well-known 
for  their  high  flexibility,  high  tensile  strength,  low  weight,  high  temperature  tolerance 
and  low  thermal  expansion.  They  are  widely  used  in  aerospace,  civil  engineering, 
military,  and  automobile.  For  Region  2  we  select  lightweight,  flame  resistant  Nomex 
honeycomb  as  the  core  material.  The  corresponding  thermal  properties  of  carbon  fibers 
and  Nomex  honeycomb  are  given  in  Table  1  (Wang  et  al.  2007;  Silva  et  al.  2005). 

It  is  worthwhile  to  note  that  the  skin  materials  have  much  larger  thermal 


conductivity  than  the  core  materials. 


582 


HONG  ZHOU,  JOSIAH  PICKETT,  AND  AURELIA  MINUT 


Table  1:  Thermal  properties  for  carbon  fibers  and  Nomex  honeycomb 


Materials 

Heat  Capacity 

(or  Specific 

Heat)  Cp 

(unit:  J/(kg*K)) 

Density  p 

(unit:  Kg  i  m3 ) 

Thermal 

Conductivity 

Kt 

(unit: 

W/(m*K)) 

Thermal  Diffusivity 

Kt 

aT=P-c„ 

(unit: m2  / s) 

Carbon 

fibers 

676 

1760 

84 

7.0602  xl0“5 

Honeycomb 

1172  ~  1340 

1100 

0.175 

1.1872xl0~7  ~1.3574xl0”7 

3.  NUMERICAL  RESULTS 

In  this  section  we  investigate  the  effects  of  different  design  parameters  on  the 
maximum  temperature  rise  of  a  sandwich-structured  composite.  We  apply  the 
commercial  software  COMSOL  to  solve  the  partial  differential  equations  (l)-(4)  with 
the  prescribed  initial  and  boundary  conditions.  COMSOL  is  a  multi-physics 
commercial  software  based  on  finite  element  methods  which  can  be  used  to  solve  our 
model  equations  efficiently. 

We  first  construct  a  sandwich  structure  where  the  core  has  thickness  d2  (mm) 
and  cross  section  dimension  304.8mm  x  304. 8mm  and  the  top  and  bottom  skins  are 
with  thickness  dx  =  di  (mm).  In  our  study  we  will  fix  d2  and  vary  di  and  then  we 

will  fix  <fj  and  vary  d2 . 

We  start  with  a  stationary  Gaussian  beam  applied  on  the  top  surface  of  the  sandwich 
structure  where  70  =500W,d  =  0.0025m,  a  =  b  =  0,  Lx  =  Ly  =  0.3048m  in  (3).  In 

Figure  2(a)  we  take  a  snapshot  of  the  temperature  rise  distribution  of  the  composite 
structure  after  the  laser  beam  is  shone  for  6  seconds  on  the  top  surface  whereas  in 
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Figure  2(b)  we  show  the  temperature  rise  in  a  vertical  slice  across  the  center  of  the  top 
surface. 


(a)  (b) 

Figure  2:  The  temperature  rise  of  a  sandwich  structure  induced  by  a  stationary 
Gaussian  beam,  (a)  a  3-D  view;  (b)  a  vertical  slice  through  the  center  of  the  top 
surface. 

From  Figure  2  we  see  that  the  maximum  temperature  (around  3358K)  occurs  at  the 
center  of  the  top  surface  which  is  exactly  the  spot  where  the  laser  beam  hits  the 
structure.  In  Figure  3  we  plot  the  transient  temperature  rise  at  the  center  of  the  top 
surface.  As  time  goes  by,  the  spot  gets  hotter  and  hotter.  Here  the  thickness  of  the  core 

is  fixed  as  d2  =  6mm  and  the  skin  thickness  is  chosen  as  dl=d3=  0.5  mm. 


s«'«tvc  •(  •  ps*«ii 


Figure  3:  The  transient  temperature  rise  of  a  sandwich  structure  at  the  center  of 
the  top  surface. 
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Now  we  vary  the  thickness  of  the  skin  from  Omm  (no  skin)  to  1mm  and  compute 
the  maximum  temperature  rise  in  the  composite  structure.  The  maximum  temperature 
rise  versus  the  thickness  of  the  skin  is  plotted  in  Figure  4(a). 


12000 

10000 

8000 
1  6000 
4000 

2000 

°0  02  0.4  0.6  0.8  1 

thickness  of  the  skin  x  -j  q  3 

(b) 

Figure  4:  (a)  The  maximum  temperature  rise  of  a  sandwich  structure  as  a 
function  of  the  skin  thickness,  (b)  The  data  in  Figure  4(a)  fitted  by  two  different 
least  squares  functions. 


(a) 


Note  that  even  with  a  thin  layer  of  skin  on  the  top  and  bottom  of  the  core,  the 
maximum  temperature  rise  can  be  reduced  significantly  in  comparison  with  the  case 
where  there  is  no  skin  at  all.  So  even  with  a  thin  layer  of  skin  one  gains  big  protection. 
In  Figure  4(b)  we  fit  the  data  in  Figure  4(a)  by  two  different  sets  of  least  squares  fitting 
functions.  The  dashed  curve  corresponds  to  the  fitting  function 

TmJhled  =  1^1  +  1233.3  whereas  the  solid  curve  is  described  by 
dx 

TmJ'“‘ed  _  0-000047  +  T5  +  q  ^he  maximum  temperature  rise  is  fitted  better  by 
dx  dx 

the  second  fitting  function  which  is  a  quadratic  function  of  the  reciprocal  of  the  skin 
thickness. 

A  practical  interpretation  of  Figure  4  can  be  as  follows.  If  we  want  to  build  an 
object  with  sandwich  structured  composite  materials  which  can  survive  laser  attacks 
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(with  power  500W)  for  six  seconds  with  a  maximum  temperature  rise  less  than,  say, 
4000K,  then  the  skin  thickness  should  be  greater  than  0.4mm. 

Another  possible  way  to  enhance  damage  tolerance  is  to  increase  the  core  thickness 
while  fixing  the  skin  thickness.  Intuitively,  increasing  core  thickness  will  enlarge  the 
bulk  volume  of  the  materials  and  therefore  is  expected  to  reduce  the  maximum 
temperature  rise  at  a  fixed  time.  We  carried  out  some  preliminary  simulations  and  our 
results  indicate  that  changing  the  core  thickness  is  not  significant  in  reducing  the 
maximum  temperature  rise. 

Now  we  would  like  to  see  how  the  beam  size  affects  the  maximum  temperature  rise 
of  the  composite.  We  fix  the  composite  structure  so  that  the  top  and  bottom  skin  layers 

have  thickness  dl=d3  =  0.5  mm  and  the  middle  core  has  thickness  d2  =6  mm  while 

varying  the  effective  beam  radius  d  from  1mm  to  5mm.  Figure  5(a)  depicts  the 
maximum  temperature  rise  as  a  function  of  the  effective  beam  radius  whereas  Figure 
5(b)  shows  two  sets  of  least  squares  fittings  with  the  data  in  Figure  5(a). 


(a)  (b) 

Figure  5:  (a)  The  maximum  temperature  rise  of  a  sandwich  structure  versus  the 
effective  beam  radius  of  a  stationary  Gaussian  beam,  (b)  The  data  in  Figure  6(a) 
fitted  by  two  different  least  squares  functions. 
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More  specifically,  in  Figure  5(b)  the  dotted  curve  is  fitted  by  the  function 
4  1 

TmJ'“ed  =  — 1 — h  1637.3  while  the  solid  curve  is  fitted  by  the  function 
d 

_  fltted  -0.0064  9.5590  , 

TmJ  = - - - h - +  646. 1341 .  It  appears  that  the  solid  curve  fits  the  data 

d  d 

better  than  the  dotted  curve.  Figure  5  shows  that  the  maximum  temperature  rise 
decreases  when  the  beam  size  increases.  This  is  expected  since  the  beam  becomes  less 
focused  when  its  size  increases. 

Next  we  study  the  temperature  rise  induced  by  a  dithering  laser  beam.  In  formula  (3) 
we  choose  I0=500W,d  =  0.0025m,a  =  0,b  =  0.05m,Lx=Ly=0.3048m,T  =  ls  .  In 

Figure  6  we  give  a  snapshot  of  the  temperature  rise  at  final  time=6s  of  the  top  surface 
(left  panel)  and  the  vertical  slice  passing  through  the  center  of  the  top  surface  (right 
panel).  Due  to  the  dithering  nature  of  the  beam,  the  maximum  temperature  rise  occurs  at 
different  spots  at  different  times.  The  maximum  temperature  rise  induced  by  a  dithering 
beam  is  dropped  significantly  to  862K  compared  with  the  stationary  beam  case  where 
the  maximum  temperature  rise  is  around  335 8K. 


(a)  (b) 

Figure  6:  The  maximum  temperature  rise  of  a  sandwich  structure  induced  by  a 
dithering  laser  beam.  Left  panel:  top  surface;  right  panel:  a  vertical  slice. 

We  also  follow  the  transient  temperature  rise  at  a  fixed  point  (the  center  of  the 
top  surface).  The  result  is  shown  in  Figure  7.  Due  to  the  fact  that  the  beam  moves  back 
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and  forth  along  a  line  and  passes  through  the  center  periodically,  the  temperature  rise  at 
a  fixed  point  oscillates  and  the  envelope  of  the  amplitude  increases  with  time. 


Figure  7:  The  transient  dithering-laser-beam-induced  temperature  rise  of  a 
sandwich  structure  at  the  center  of  the  top  surface. 

For  a  dithering  laser  beam,  we  will  investigate  the  effect  of  the  skin  thickness  as 
well  as  the  effect  of  the  beam  size  on  the  maximum  temperature  rise.  From  the  previous 
study  of  the  stationary  beam,  the  maximum  temperature  rise  is  not  very  sensitive  to  the 
core  thickness  and  therefore  will  not  be  studied  any  more.  In  Figure  8(a)  we  show  the 
maximum  temperature  rise  dependence  on  the  thickness  of  the  skin  layer  while  in 
Figure  8(b)  we  fit  the  data  with  two  different  least  squares  functions.  The  dotted  curve 

0  2210 

corresponds  to  the  fitted  function  TmsJltted  =— - +  401.7029  whereas  the  solid 

dx 

■a  u  a  u  t1  nued  -0.000016,0.395,^^  A 
curve  is  described  by  TmsaJ  = - - - h - 1-143.20  .  As  expected,  a 

d}  dx 

combination  of  sandwich  structure  and  dithering  beam  leads  to  a  big  reduction  in  the 
maximum  temperature  rise. 
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(a)  (b) 

Figure  8:  (a)  The  dithering-beam-induced  maximum  temperature  rise  of  a 
sandwich  structure  as  a  function  of  the  skin  thickness,  (b)  The  data  fitted  by  two 
different  least  squares  functions. 

In  Figure  9(a)  we  plot  the  maximum  temperature  rise  as  a  function  of  the 
effective  laser  beam  radius.  As  the  beam  size  increases,  the  maximum  temperature  rise 
decreases.  We  fit  the  data  in  Figure  9(a)  by  two  different  functions  in  Figure  9(b).  Here 
the  dashed  curve  is  described  by  a  linear  function  of  the  reciprocal  of  the  beam  radius 

Tjitted  =  +  56  j  .4320  and  the  solid  curve  is  given  by  a  quadratic  function  of  the 

d 

reciprocal  of  the  beam  radius  TmJ,tted  =  ^'^^  +  ^'^^^  +  275.2048.  Clearly,  the 

d  d 

maximum  temperature  rise  behaves  more  like  a  quadratic  function  of  the  reciprocal  of 
the  beam  size. 


(a)  (b) 

Figure  9:  (a)  The  maximum  temperature  rise  of  a  sandwich  structure  versus  the 
effective  beam  radius  of  a  dithering  Gaussian  beam,  (b)  The  data  fitted  by  two 
different  least  squares  functions. 
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Finally  we  move  to  the  temperature  rise  induced  by  a  rotating  laser  beam.  First 
in  Figure  10(a)  we  take  a  snapshot  of  the  temperature  rise  induced  by  a  rotating  laser 
beam  with  I0  =  500W,  d  =  0.0025m,  a  =  b  =  0.05 m,Lx  =  L  =  0.3048m  in  expression 

(3). 


(a)  (b) 

Figure  10:  The  maximum  temperature  rise  of  a  sandwich  structure  induced  by  a 
rotating  laser  beam.  Left  panel:  top  surface;  right  panel:  a  vertical  slice. 

Figure  10(a)  indicates  that  the  maximum  temperature  rise  drops  further  to  about 
323K.  Figure  10(b)  shows  the  temperature  rise  on  a  vertical  slice. 

In  Figure  1 1  we  plot  the  transient  temperature  rise  at  the  center  point  of  the  top 
surface.  Here  the  temperature  rise  increases  monotonically  as  a  function  of  time  even 
though  the  magnitude  of  the  temperature  rise  is  very  small.  The  increasing  nature  of  the 
temperature  is  due  to  the  insulating  boundary  condition. 
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transient  temperature  at  a  fixed  point 


Figure  11:  The  transient  rotating-laser-beam-induced  temperature  rise  of  a 
sandwich  structure  at  the  center  point  of  the  top  surface. 

Now  we  want  to  see  how  the  maximum  temperature  rise  behaves  as  we  vary 
either  the  skin  thickness  or  the  effective  beam  size.  Our  numerical  results  are  shown  in 
Figures  12  and  13,  respectively,  together  with  least  squares  fitting  functions. 


(a)  (b) 

Figure  12:  (a)  The  rotating-beam-induced  maximum  temperature  rise  of  a 
sandwich  structure  as  a  function  of  the  skin  thickness  dx .  (b)  The  data  in  Figure 

12(a)  fitted  by  two  different  least  squares  functions. 
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0  0469 

In  Figure  12,  the  dashed  curve  is  given  by  Tm,Jitted  =  — - 1-210.4757 

dx 

,  .  r,  .  „  fMed  -0.000005  0.104  __ 

whereas  the  solid  curve  is  T  j  = - - - -i - +  125.98 . 

max  7  z  7 

dx  d] 

Figure  12  and  Figure  13  again  demonstrate  that  the  maximum  temperature  rise  can 
be  described  by  a  quadratic  function  of  the  reciprocal  of  the  skin  thickness  or  the  beam 
size  where  the  coefficients  in  the  quadratic  form  are  different  for  each  case. 


500i - . - . - . - . - . - , - , - 1  420 


(a)  (b) 


Figure  13:  (a)  The  maximum  temperature  rise  of  a  sandwich  structure  versus  the 
effective  beam  radius  of  a  rotating  Gaussian  beam  d .  (2)  The  data  in  Figure  13(a) 
fitted  by  two  different  least  squares  functions.  The  dashed  curve  is  described  by 


fitted  _ 


0.35 


+  173.65  and  the  solid  curve  is  T 


fitted  _ 


-0.000527  0.7998 

- - - + - +  92.01. 


4.  CONCLUDING  REMARKS 


We  have  shown  numerically  how  the  maximum  temperature  rise  induced  by  a 
laser  beam  for  a  composite  sandwich  structure  scales  as  a  function  of  the  thickness  of 
the  skin  layer  or  the  effective  beam  size.  Our  study  includes  three  kinds  of  laser 
beams:  a  stationary  Gaussian  beam,  a  dithering  Gaussian  beam  or  a  rotating  Gaussian 
beam.  All  the  cases  share  a  universal  rule:  the  maximum  temperature  rise  behaves  as  a 
quadratic  function  of  the  reciprocal  of  the  skin  thickness  or  the  effective  beam  size. 
Validation  of  our  predictions  with  the  experimental  data  will  be  pursued  in  the  future. 
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